Abstract. We present the first optical images of scat- 
tered light from large, detached gas/dust shells around 
two carbon stars, R Scl and U Ant, obtained in narrow 
band filters centred on the resonance lines of neutral K 
and Na, and in a Stromgren b filter (only U Ant). They 
confirm results obtained in CO radio line observations, but 
also reveal new and interesting structures. Towards R Scl 
the scattering appears optically thick in both the K and 
Na filters, and both images outline almost perfectly circu- 
lar disks with essentially uniform intensity out to a sharp 
outer radius of «21". These disks are larger - by about 
a factor of two - than the radius of the detached shell 
which has been marginally resolved in CO radio line data. 
In U Ant the scattering in the K filter appears to be, at 
least partially, optically thin, and the image is consistent 
with scattering in a geometrically thin («3") shell (radius 
Ri43") with an overall spherical symmetry. The size of this 
shell agrees very well with that of the detached shell seen 
in CO radio line emission. The scattering in the Na filter 
appears more optically thick, and the image suggests the 
presence of at least one, possibly two, shells inside the 43" 
shell. There is no evidence for such a multiple-shell struc- 
ture in the CO data, but this can be due to considerably 
lower masses for these inner shells. Weak scattering ap- 
pears also in a shell which is located outside the 43" shell. 
The present data do not allow us to conclusively identify 
I the scattering agent, but we argue that most of the emis- 
O ^. sion in the K and Na filter images is to due to resonance 
line scattering, and that there is also a weaker contribu- 
tion from dust scattering in the U Ant data. Awaiting new 
■ observational data, our interpretation must be regarded as 
Ci • tentative. 
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1. Introduction 

Low, and intermediate-mass stars populate the asymp- 
totic giant branch (AGB) after central He-burning. During 
this stage of evolution they are known to experience in- 
tense mass loss in the form of a cool, slow, and essentially 
spherically symmetric wind (Olofsson 1996 ). The mass loss 
eventually becomes so large that it, rather than the nu- 
clear burning, determines the evolutionary time scale of 
these stars. This limits the AGB life time, and e.g. pro- 
vides an explanation to the discrepancy between the ob- 
served and calculated luminosity functions of stars on the 
horizontal branch and on the AGB in globular clusters 
Ibcn 



(Rcnzini 1981 



1981). The mass loss also makes an 



important contribution to the enrichment of the interstel- 
lar medium in terms of nu clear- processed material and 
dust particles (Busso et al. 1999). Considering this, it is 



unfortunate that our knowledge of the details of the mass 
loss and the mechanism behind it is so limited. 

The formation of circumstellar envelopes (CSEs) of gas 
and dust around AGB stars is a consequence of the mass 
loss. The main constituent of the outflowing gas, H2, is 
generally not observable, despite being excited by shocks 
in a few particular circumstances, notably young post- 
AGB objects an d plan etary nebulae (Beckwith et al. 1978 ; 
Weintraub et al. 1998 ). Standard stellar atmosphere chem- 
istry theory predicts that the first molecule to form as the 
gas cools is CO (Glassgold 1996 ). Depending on the star 
being 0-rich (i.e., [C]/[0]<1) or C-rich ([C]/[0]>1) the 
CO molecules consume all the available carbon or oxygen 
atoms, respectively. This leads to a drastically different 
molecular chemistry and dust composition in the CSEs 
around these two types of AGB stars. Direct measures of 
the gas, the principal component in terms of mass of a 
CSE, are provided by radio line emission of these minor 
molecular constituents. Of particular importance for the 
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study of the mass loss properties are observations of CO 
for both O- and C-rich stars (Olofsson et al. 1996| ). 

The temporal variation of the mass loss rate is to a 
large extent unknown. This applies to all time scales from 
the pulsation period to the full time scale for the AGB- 
phase. For the shortest time scales we are limited by the 
spatial resolution of the observations, while for the longest 
time scales we lack suitable observational probes. On the 
intermediate time scales (10^-10"* yr) there is now grow- 
ing evidence for substantial variations in the mass loss 
rate, e.g., detached CO and dust shells have been detected 
around a number of AGB and post-AGB stars (Olofsson 
et al |1990| , |1996| ; Lind gvist et al. |1996| |1999| ; Waters et ah 



1994; Izumiura et al. 1996, 1997; Hashimoto et al. 1998 



Speck et al. 20001 ), and multiple-shell structures have been 



seen in scattered light towards IRC+10216 and a num- 



ber of post-AGB objects (Harpaz et al. 1997; Kwok et al 



1998; Sahai et al. 1998; Mauron & Huggins 1999). Interfer 



ometric observations of the CO shell around the carbon 
star TTCyg provide a splendid result (Olofsson et al. 
1998i I2OOOI) . The sheU is large (radius w35"), geometri- 
cally thin (average width «2.5"), and remarkably spher- 
ical (only ±3% variation in radius). Similar results have 
been obtained for the carbon star U Cam (Lindqvist et al. 
1999 ). A significant fraction of the AGB-stars has excess 



emission at 60 /im (and 100/im) in the IRAS PSC, which 
has been interpreted as a sign of detached envelopes, and 
hence va riable mass loss ( Willem s & de Jong 1988; Chan 
& Kwok 199C; Zijlstra et al. 1992). However, in most cases 



this excess can probably be attributed to interstellar cir- 
rus emission (Ivezic & Elitzur 1995). Thus, there remain 
only a few, clear examples of highly episodic mass loss 
on the AGB. The exact origin of the detached shells is 
still uncertain. They could be rare examples of 'thermal 
pulse '-induced mass loss eruptions (Olofsson et al. 1990; 
Schroder et al. 1999), but effects of interacting winds can- 



not be excluded, with the CO shells being the shock zones 
(Olofsson et al. |2000t SteflFen & Schonberner |200C| ). 

R Scl and U Ant are among the few carbon stars with 
extended, detached shells detected in CO radio lines (Olof- 
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sson et al. 1990, 1996). For R Scl, the CO brightness maps 
show a marginally resolved detached shell (radius ~10"), 
which appears to have detached fairly recently (^10'^ yr). 
Detections of other molecules, such as HCN, CN and CS, 
point to a dense shell where these molecules are still rel- 
atively abundant and excited (Bujarrabal & Cernicharo 
1994| ; Olofsson et al. |1996D . The CO emission from U Ant 
shows both a 'normal' CSE and a detached shell (radius 
«40"). The latter appears to have an overall spherical 
symmetry, but its brightness distribution is not smooth, 
suggesting dumpiness at some level in the density distri- 
bution, as well as a marked decrease in the CO inten- 



sity towards the south-west. Izumiura et al. ( 1997 ) decon- 
volved IRAS 60 and 100 /im observations of U Ant, and 
they inferred the existence of two extended dust shells. 
The (uncertain) size of the inner shell, apparent as a 
brightness distribution broader than the PSF, is consis- 
tent with that of the detached CO shell. The outer shell 
has an inner radius of about 3' and no CO counterpart. 

The interpretation of the CO emission is hampered by 
the fact that it depends on the excitation as well as on the 
chemistry (particularly the photodissociation by the in- 
terstellar radiation field) of the CO molecules. Hence, the 
relationship between the density distribution (and conse- 
quently the mass loss history) and the CO brightness dis- 
tribution is uncertain. Also, the lack of an interferometer 
in the southern hemisphere prevents us from obtaining CO 
radio line observations of R Scl and U Ant with the high 
resolution achieved for the northern sources TT Cyg (Olof- 
sson et al. |1998| , pOOOD and U Cam (Lindqvist et al. |1996| , 
1999). Infrared emission from dust is the only probe of 



CSEs which is not destroyed by the interstellar radiation 
field and therefore these observations trace the long-term 
history of mass loss. However, the data are often of very 
poor spatial resolution. Hence, there is an urgent need for 
other ways of observing the detached shells. In this pa- 
per we present the first optical images of the envelopes 
around R Scl and U Ant, which we infer to be the result 
of stellar light scattered in atomic resonance lines and/or 
by dust particles. These high spatial resolution data will 
allow us to further study the morphology and structure of 
the shells. The feasibility of optically observing CS Es wa s 
first shown spectroscopically by Bernat & Lambert ( |1975|) . 
They used a method whereby normalized on-star spectra 
are subtracted from off-star spectra, to reveal the faint 
circumstellar emission; this technique was used and devel- 
oped further by Mauron & Caux (1992), Plez & Lambert 
(|1994D, and Gustafsson et al. (|1997D. 



2. Observations and data reduction 

2.1. Observations 

The imaging observations presented in this paper were 
done primarily in narrow (5 nm FWHM) filters centred 
on the resonance lines of Nal (589.0 nm and 589.6 nm; 
the D fines) and KI (769.9 nm), using the EFOSCl and 



Table 1. Observations performed during the two observ- 
ing runs [ND, no shell emission detected; D(l), shell emis- 
sion detected during the first (test) run; D(2), shell emis- 
sion detected during the second run] 



KI 



NaD Str.b Str.v Gunnz 



R Scl 


D(l) 


D(l) 




U Ant 


D(2) 


D(l,2) D(2) ND 


ND 


S Set 


ND 


ND 




U Hya 


ND 


ND 




R Hya 


ND 


ND 




X TrA 


ND 


ND 




VX Sgr 


ND 


ND 





EF0SC2 focal reducer cameras on the ESO 3.6m telescope 
in 1994 and 1999, respectively. We used narrow filters to 
reduce the amount of direct photospheric light in the im- 
ages, and thereby increase the contrast of any circumstel- 
lar line scattered light. Even narrower filters were consid- 
ered, but they would introduce new unwanted difficulties 
in the observations, such as the temperature dependence 
of the transmitted wavelengths, and wavelength shifts over 
the image that are greater than the filter widths (pro- 
duced by the large difference in angle at which light passes 
through the filter from one side of the field to the other in 
a parallel beam instrument). During the second observing 
run, we also observed U Ant in other filters, which con- 
tain no strong resonance lines (Stromgren v, A=410.0nm, 
AAf=18nm; Stromgren b, A=468.0nm, AAt=20nm; and 
Gunnz, A=900nm, AAf=20nm, approximately). In Ta- 
ble |l| we summarize the observations done during the two 
observing runs. 

The stellar fluxes turned out to be typically four orders 
of magnitude higher than the circumstellar scattered light 
fluxes. A coronograph is therefore essential, and we used 
an opaque dot placed on a glass plate in the focal plane 
of the EFOSC instruments. By using this method, most 
of the photospheric light that is scattered in the Earth's 
atmosphere is avoided. However, even the use of a corono- 
graphic mask is not enough, and an additional method is 
required in order to subtract the remaining direct photo- 
spheric radiation. This additional method is based on the 
observation of template stars, i.e., stars with about the 
same magnitude and spectral type as the target stars, but 
without CSEs. These template star images are subtracted 
from the target star images during the reduction. Due to 
atmospheric variations, the stellar profiles change during 
the night. Therefore, we observed several template stars 
well distributed throughout the full runs. Standard stars 
were only observed during the second observing run. 

The first images in circumstellar scattered light were 
obtained during observations carried out in 1994 at a test 
run with the ESO 3.6m telescope and the EFOSCl 'fo- 
cal reducer'-type camera. The pixel plate scale was 0'.'61, 
and exposure times of only a few minutes were sufficient 
for the detections. Three bright carbon stars with de- 
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tached shells detected in CO radio lines (R Scl, U Ant, 
S Set) were chosen as targets during this test run. Emis- 
sion due to circumstellar scattering in extended envelopes 
was found for two of them, R Scl and U Ant. These de- 
tections are, however, of relatively low S/N-ratio. Around 
S Set we found, using a polarimetric technique, an ex- 
tremely faint and diffuse shell at very low S/N-ratio. This 
test run turned out to be the first success after a series of 
attempts in previous years with other telescopes [the 2.5m 
Nordic Optical Telescope (NOT) and the MPG/ESO 2.2m 
telescope] . Successful spectroscopic observations were per- 
forr aed w ith the ESO CAT/CES-system by Gustafsson et 
al. ( 1997 ) demonstrating the presence of circumstellar res- 
onance line scattering around R Scl and other stars. The 
observational problem of imaging the circumstellar shells 
is clearly their faintness, which makes the stellar light scat- 
tered in the Earth's atmosphere and the telescope very 
problematic. The ESO 3.6m telescope was crucial to the 
outcome of the project. Apart from a large collecting area, 
it has an equatorial mounting. In alt-az telescopes, like the 
NOT and the ESO NTT, the rotation of the instrument 
during the observations smears out the spider diffraction 
pattern over the images. This results in a substantial loss 
of azimuthal information, rendering it difficult to make a 
detection of weak emission from an extended envelope. 

Based on the successful observations during the test 
run, we continued in 1999 with the imaging of one of the 
stars, U Ant, with the EF0SC2 camera. The pixel plate 
scale in EF0SC2 is 0'.'32. We detected again the emission 
from the envelope around U Ant, in this case in several fil- 
ters and with much higher S/N-ratio. During this second 
run, we also observed two stars with detached dust shells, 



U Hya (Waters et al., |199J) and R Hya [the only 0-rich 
sou rce w ith a detected detached dust shell (Hashimoto et 
al., 1998| )], and one source wit h a p ossible detached dust 
shell, X TrA (Izumiura et al. 1995 ). In addition, we ob- 
served VX Sgr, an M-type supergiant. We were not able 
to detect any circumstellar emission in these sources. Due 
to the complexity of the observations, there can be sev- 
eral reasons for this negative result. Normal envelopes, in 
which the density peaks close to the star, are very difficult 
to detect in scattered stellar light. It is worth noting that 
the circumstellar emission observed towards IRC-l-10216 
is due to scattered interstellar light, and the central star 



is highly obscured (Mauron & Huggins 1999). Shells lo- 
cated close to the star, more extended (and therefore even 
fainter) shells (e.g., in S Set), and the presence of the 
star in a crowded stellar field are other plausible reasons 
for the non-detections. We postpone the discussion of the 
tentative S Set shell and the non-detections to a future 
paper. 

2.2. Data reduction 

The image reduction was done in two steps using IRAF. 
First, an ordinary CCD reduction (bias subtraction, flat- 





Distance from star centre (arcsec) 



Distance tram star centre (arcsec) 



Fig. 1. Azimuthally averaged radial brightness profiles of 
a target star (full fine) and a template star (dot-dashed 
line) in the KI filter. In the case of a shell detection (U Ant, 
left), the target star profile has an additional component 
from the light scattered in the shell. In the case of no 
detection (HD184835, right), the target and template star 
radial profiles match well 

field correction, and cosmic ray removal) was performed. 
This resulted in images with some considerable distor- 
tions. In particular, the coronographic optics at the ESO 
3.6m telescope was fairly crude, e.g., the focal plane ob- 
scuring spot masks are made on a glass plate that in- 
troduced optical distortions. Furthermore, the lack of an 
optimized apodising mask to remove the spider diffrac- 
tion pattern on the images, resulted in two highly satu- 
rated stripes, corresponding to the segments of the spider, 
along the horizontal and vertical directions of the images, 
where the reduction is much more complicated. Therefore, 
the circumstellar emission detected in the various images 
shown in this paper is not reliable along these two di- 
rections. Some cleaning was performed in the images in 
order to remove field stars whose brightness could swamp 
the circumstellar emission. 

Second, the template star subtraction was done (the 
template star images are first reduced as above and then 
median filtered in this process). One of the key aspects 
of these observations is the placement of the objects un- 
der the coronographic spot. This is a very difficult task 
and when doing the reductions we have realized that the 
stars were not positioned exactly under the centre of the 
coronographic mask, and therefore the positions of the 
stars are difficult to locate precisely. An accurate subtrac- 
tion requires the target and template stars to be perfectly 
aligned, and so an indirect method had to be used to 
achieve this. We estimated the positions of the stars us- 
ing their radial profiles along various position angles. We 
defined these positions as the middle points between two 
pixels (one on each side of the centre) with the same inten- 
sity level. We averaged the values obtained using different 
intensity levels and different radial profiles (along different 
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position angles) to get a final location of the target and 
template stars. 

Next in the reduction came the normalization of the 
template star images, i.e., normalizing the template star 
intensity to that of the target star. Again, the absence of 
peaks in the images made this task somewhat complicated. 
We normalized by a method of trial and error, using a ra- 
dial range outside the innermost usable points and 
5" from the star centres in the reductions of the R Scl and 
U Ant images, respectively) of several radial profiles. The 
profiles of the target star and the template star are usually 
perfectly matched close to the coronographic mask after 
normalization. The addition of a small constant value to 
the template image was used to equalize the background 
light in the images due to the sky background. The in- 
ner points, where the intensity is much higher, are not 
significantly affected by this addition. 

Finally, we subtracted the template star image from 
the target star image. This subtraction removes both the 
stellar and the background components from the target 
star image, leaving us only with the circumstellar emis- 
sion. We have repeated this procedure using two template 
stars, in order to make sure that the observed circumstel- 
lar contributions are not artifacts produced in the tele- 
scope or during the reduction. As expected, in this case 
only noise is left. Fig. |l| shows two examples of azimuthally 
averaged radial profiles in the KI filter: a circumstellar 
shell detection and a non-detection. In the first case, the 
radial profile of the star (U Ant) clearly shows additional 
emission when compared to the template star profile. In 
the case of a non-detection, both profiles are essentially 
the same, and their subtraction produces only noise. 

Due to the fact that the centring and normalization 
cannot be performed exactly, the high values, combined 
with the large gradient, close to the star make it impos- 
sible to follow any shell emission inside a certain radius, 
and we have indicated this area in the final images using 
a central, black spot. In the following we will only discuss 
the properties of the scattered emission outside this sub- 
jectively chosen area. We also note here that the template 
star subtracted (TSS) images and radial profiles become 
progressively less reliable the closer the angular offset gets 
to this area. 

It was not possible to calibrate the R Scl images in 
absolute fluxes because no standard stars were observed 
during the test run. However, we have done an estimate 
of the shell fluxes by considering all possible effects which 
determine the detected number of ADUs in the CCD, see 
Sect. 3.4. In the case of U Ant, we have done a proper 
calibration of the shell fluxes in the resonance filters, us- 
ing observations of a tertiary spectrophotometric standard 
star (LTT 3864) taken from the sample by Hamuy et al. 

O)- 

3. Results 
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Fig. 3. Azimuthally-averaged radial brightness profiles of 
the circumstellar emission around R Scl in the KI (left) 
and NaD (right) filters. The inner parts of the profiles are 
affected by the presence of the coronographic mask and 
the template star subtraction and are therefore not shown 
here. The dashed line shows the best-fit seeing-smoothed 
step function (see text for details) 

3.1. TSS images of R Scl in the KI and NaD filters 

The TSS images of R Scl in the KI and NaD filters re- 
veal remaining brightness distributions in the form of es- 
sentially uniform-intensity disks outside the masked ar- 
eas. Fig H The NW regions of the disks show a some- 
what enhanced emission in both filters, but there remains 
distortions in the images that limit the useful azimuthal 
information. The azimuthally averaged radial brightness 
profiles are relatively constant (outside the masked area) 
within an outer radius of «21" in both filters (correspond- 
ing to 1.1x10^ ''cm at a distance of 360 pc, see Sect. 5.2), 
Fig. ||. The outer radius is obtained by fitting a step func- 
tion, smoothed by a seeing Gaussian, to the radial pro- 
files, (see Sect. 3^). The outer radius is defined as the 
half power radius of the fit. The outer edges of the ra- 
dial brightnesses are not as sharp as the seeing smoothed 
step function, but the half power radii are only larger by 
w2" than the radii where the intensities start to drop. 
The emission is fainter in the NaD filter. We interpret the 
relative faintness in the Na D image as due to the star be- 
ing fainter in this filter, and the scattering being at least 
partially optically thick (see Sect. ^). 

In Fig. ^ we plot the estimated disk radii (in the KI 
and Na D filters) , defined again as the half power radius of 
a fitted smoothed step function, as a function of position 
angle (PA). The brightness disks appear close to circular 
in both filters, with only small deviations at the ^5% level. 
This indicates an emitting region with an overall spherical 
symmetry. The deviations can be partly attributed to the 
uncertainty in the position of the star behind the corono- 
graphic mask, but there is no clear tendency for an overes- 
timate of the radii in a range of PAs and an underestimate 
in the complementary interval, as would be expected from 
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Fig. 2. Template star subtracted images of R Scl in the KI (left) and NaD (right) filters obtained with EFOSCl. The 
black disks mask the regions where reliable data on the scattered light does not exist 



onance line filters (Olofsson et al. 1996). We will further 
discuss this in Sect. t).2l 



1 OO 200 
Position ongio 



Fig. 4. Outer radii (defined in the text) of the R Scl cir- 
cumstellar emission in the KI (upper) and NaD (lower) 
filters measured at different position angles. The thick line 
gives the mean radius obtained from these points. The er- 
ror bar indicates the uncertainty in the radius estimates 
due to the pixel size 



a bad centring. Furthermore, the deviations seem to be re- 
produced in the same manner in both filters, suggesting 
actual slight departures from spherical symmetry. 

Deconvolution of CO(J=3-2) data, obtained with a 
16" beam, using the Maximum Entropy Method shows a 
marginally resolved detached shell with a CO peak inten- 
sity radius which is significantly smaller, by about a factor 
of two, than the scattering envelope observed in the res- 



3.2. TSS images of U Ant in the KI and NaD filters 

The TSS images of U Ant in the KI and NaD filters pro- 
vide spectacular results. Fig ^. In the KI filter image a 
geometrically thin, almost circular, ring with a diameter 
of «85" is centred on the star. This is most likely the re- 
sult of limb brightening in a geometrically thin, essentially 
spherical shell. Due to the bad coronographic optics, the 
image contains striation and diffraction spikes that could 
not be removed during the reduction, and this limits the 
available azimuthal information. 

The azimuthally-averaged radial brightness profile in 
the KI filter presented in Fig. ^ shows a very pronounced 
intensity peak, and the mean shell radius was estimated 
to be 43" (the corresponding linear radius is 1.7x10^^ cm 
at the adopted distance of 260 pc, see Sect. ^ ) by fit- 
ting a brightness distribution appropriate for an optically 



thin homogeneous shell (see Sect. 3.4). This shell radius 
estimate from the scattered light agrees very well with the 
peak intensity radius of the detached CO shell observed by 
Olofsson et al. (1996), w41" (they have estimated their an- 
gular resolution to be wl5", but the peak intensity radius 
has an uncertainty of only a few arc seconds). The emission 
tapers off gradually towards the centre, and there are in- 
dications of additional structure, which is more prominent 
in the Na D filter data. The emission decreases drastically 
beyond the peak radius, but the radial profile also shows 
a small bump at «50". 
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Fig. 5. Template star subtracted images of U Ant in the KI (left) and NaD filters obtained with EF0SC2. The black 
disks mask the regions where reliable data on the scattered light does not exist 



In the TSS image of U Ant in the NaD filter we also 
clearly see the circumstellar emission, Fig. ^. However, the 
morphology is quite different from that obtained in the 
KI filter, in the sense that the brightness stays at nearly 
constant level much closer to the star. In fact, the ap- 
pearance suggests a multiple-shell structure, which is also 
seen in the azimuthally averaged radial brightness profile. 
Fig. 1^ At least two, maybe three, different components 
are discernible, and the two outer components may not be 
concentric. Fitting sh ell b rightness distributions to these 
components (see Sect. 3.4), we have determined shell radii 
of roughly 25", 37", and 43" (1.0, 1.4, and 1.7x10^^ cm at 
a distance of 260 pc), although the result for the innermost 
component is very uncertain. The outer component coin- 
cides with the dominant peak in the KI filter image, and 
the peak radius of the detached CO shell. The difference 
between the observations in the Na D and the KI filters can 
be explained in terms of optical depth effects (see Sect. 
The possible existence of multiple shells is surprising and 
this will be further discussed in Sect. 5.2. As in the KI 



filter image, there is a small bump in the radial brightness 
profile at r^SO" offset from the star. This emission is only 
marginally visible in the image of the shell. 

Fig. H shows the estimated shell radii (for the 43" 
shell) at different PAs in the KI and NaD filter images 
of U Ant. The deviations from circular symmetry are very 
small, ;$3%. Thus, we can also infer a remarkable overall 
spherical symmetry for the U Ant shell. The small de- 
viations from the means follow the same pattern in both 
filters, indicating minor irregularities in the shell morphol- 
ogy. The KI filter image reveals edge-effects, which appear 
in the form of double-peaks in the radial brightness pro- 
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Fig. 6. Azimuthally-averaged radial brightness profile of 
the circumstellar emission around U Ant in the KI filter. 
The inner part of the profile is affected by the presence of 
the coronographic mask and the template star subtraction 
and is therefore not shown here. The dashed line shows the 
best fit shell brightness distributions (see text for details) 



files. Fig. ||. We interpret this as due to small undulations 
in the shell structure along the line-of-sight, but a more 
complex origin cannot be excluded. 

3.3. TSS images of U Ant in other filters 

We have also made observations of U Ant in filters contain- 
ing no strong resonance lines (v, b, and Gunnz). Fig. [l^ 
shows the TSS image of U Ant in the b filter. We have 
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Fig. 7. Azimuthally averaged radial brightness profile of 
the circumstellar emission around U Ant in the Na D filter. 
The inner part of the profiles is affected by the presence of 
the coronographic mask and the template star subtraction 
and is therefore not shown here. The dashed line shows the 
best-fit shell brightness distribTitions (sec text for details) 
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Fig. 8. U Ant shell radii (for the 43" shell) measured at 
different position angles in the Kl and NaD filters. The 
solid lines give the mean values for the radii obtained from 
these points. The error bars indicate the uncertainty in the 
estimates produced by the pixel size 



binned the image by 2x2 pixels in order to increase the 

sensitivity to extended emission. We find that weak emis- 
sion is detected, especially towards the south. The pres- 
ence of a bright field star with a very long horizontal 
diffraction pattern makes the detection more problem- 
atic towards the north. The azimuthally averaged radial 
brightness profile clearly shows an emission bump extend- 
ing from about 35" to 55". We have estimated its radius 
to be W46" by fitting a Gaussian to the radial brightness 
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Fig. 9. Azimuthally averaged radial brightness profiles in 
the Kl image of U Ant in PA intervals where double-peaks 
occur. The profiles have been normalized to their peak val- 
ues and shifted vertically by 0.3, 1.2, and 1.7, respectively 
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Fig. 10. Template star subtracted image of U Ant in the 
b filter obtained with EF0SC2. The black disk and stripe 
mask the region where reliable data on the scattered light 
does not exist 

profile. This is about the same radius as those of the small 
outer bumps in the Kl and Na D radial brightness profiles 
discussed above. We report also an upper limit, obtained 
from a short exposure, to any shell emission in the Gunnz 
filter. 

3.4- The stellar and circumstellar fluxes 

In order to analyse the circumstellar emission and the cir- 
cumstellar structure it is enough to compare the relative 
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Fig. 11. Azimuthally- averaged radial brightness profile of 
the circumstellar emission around U Ant in the b filter. 
The inner part of the profile is affected by the presence of 
the coronographic mask and the template star subtraction 
and is therefore not shown here 



fluxes of the circumstellar and stellar emission. Unfortu- 
nately, the use of a coronographic mask in the observations 
makes it difficult to obtain the stellar flux from an image. 
We have therefore tried to fit an extended point spread 
function of the type found by King ( |l97lD and Piccirillo 
( |1973D to estimate the stellar flux, but the resulting val- 
ues appear too low by about two orders of magnitude. We 
have not identified the reason for this, but suspect that the 
use of a coronographic mask may be part of the problem. 
Therefore, an indirect method using a synthetic spectrum 
of R Scl had to be used to estimate the stellar fluxes at 
the central wavelengths of the various filters. We adopted 
mv=5.8 and 5.5 for R Scl and U Ant, respectively. The 
variability of R Scl and U Ant (rsO.5 mag at visual wave- 
lengths) introduces some uncertainty in the calculation, 
but the obtained values should be within a factor of two 
of the actual fluxes. The results are given in Tables || and 
i 

The U Ant shell fluxes (except in Gunnz) were cali- 
brated using photometric standard stars. Sect. 2.2. In the 



case of R Scl (and U Ant in Gunn z) calibration data are 
lacking, and we have therefore converted from counts/s in 
the detector to fluxes above the atmosphere by correcting 
for atmospheric absorption, telescope aperture with a cen- 
tral obstruction, telescope reflectivity, instrument trans- 
mission, filter bandpasses and transmission, and the quan- 
tum efficiencies and gains of the CCDs. We estimate that 
the R Scl shell fluxes are accurate to within a factor of 
five, and those of U Ant within a factor of two. 

We have determined the circumstellar fluxes by fitting 
adopted brightness distributions to the azimuthally av- 
eraged radial brightness profiles, in order to be able to 
estimate the brightnesses in the occulted regions, and also 
to separate the fluxes from the different components. In 



the case of R Scl we adopt a step function, smoothed by 
a seeing Gaussian, to fit the radial profiles. The outer ra- 
dius is defined as the half power radius. Fig. ||. In the case 
of U Ant we have chosen, guided by the KI filter image, 
to assume optically thin scattering in homogeneous shells. 
This produces a generic brightness distribution with the 
shell radius and width, and the peak intensity as the fitting 
parameters. We have fitted such distributions, smoothed 
with a seeing Gaussian, to the three peaks (the compo- 
nents shelll, shell2, and shell3 in Table ||) in the NaD 
filter radial profile (ignoring for the moment the compo- 
nent outside the 43" shell). We have done the same for 
the KI filter radial profile, except that the radii of the two 
inner shells were kept the same as for the Na D filter radial 
brightness profile fit. The radial distribution of the outer 
component appears different from that of the others, and 
we have fitted a Gaussian distribution to the b filter ra- 
dial brightness profile of U Ant (the component shelly in 
Table ^). A Gaussian with the same location and width 
was subsequently fitted to the outer peaks in the KI and 
NaD data. The final fits are shown in Figs. |6[ ^ and |ll|. 

The shell surface brightnesses and fluxes obtained from 
the fits, and the ratio of circumstellar-to-stellar fiuxes 
(CS/S-ratio) are given for both stars in the different fll- 
ters in Tables || and ^. Also given are the estimated outer 
radii of the R Scl brightness distributions, and the radii 
and widths of the U Ant brightness distributions. 



4. Analytical approach 

We will in this section derive some analytical formulae 
appropriate for a simple analysis of our observational re- 
sults in terms of line or dust scattering, in the optically 
thin regime, of central stellar light in a geometrically thin, 
spherical shell. Even a rough estimate shows that in our 
cases the stellar light dominates completely over the in- 
terstellar light, as opposed to e.g. the case of IRC-l-10216 
which is heavily obscured and where the scattered light 
is of interstellar origin (Mauron & Huggins |1999| , ^OOOl ). 
Clearly, the CS/S-ratio in a given filter cannot exceed 



- AAf 

where is the ratio of the average stellar flux density 
in the total scattering wavelength range, AAgc, and the 
average stellar flux density in the fllter passband, AAf 
(e.g., aeff<l in the case of resonance line scattering where 
strong photospheric line absorption removes photons in a 
wavelength range which is much narrower than the fllter 
width). 

In the cool, detached shells of AGB stars the total line 
scattering wavelength range is determined by the turbu- 
lent motion of the gas, i.e., we assume there is no velocity 
gradient across the shell. We have adopted a local velocity 
width of lkms~^, a value which is typically assumed in 
the analy sis of CSE radio line emission (Schoier & Olofs- 
son ^001| ), and which is consistent with the CO radio line 
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Table 2. Results obtained from fitting a step function to the observed circumstellar radial brightness distributions in 
the KI and Na D filters of R Scl. The stellar fluxes are calculated (see text for details) 



Filter 


component 


Rout (") 


/ (erg/s cm^ "'') 


F (erg/s cm^) 


CS/S-ratio 


KI 


shell 


21.0 


2.0x10"'=' 


2.8x10"'^ 


8.5x10"* 




star 






3.3x10"^ 




NaD 


shell 


20.4 


1.1x10"'^ 


1.5x10"'^ 


l.lxlO-" 




star 






1.4x10"^ 





Table 3. Results obtained from fitting shell brightness distributions to the observed circumstellar radial brightness 
distributions in the KI and NaD filters of U Ant. The component external to the 43" shell has been fitted with a 
Gaussian. The stellar fluxes are calculated (see text for details) 



Filter 


component 


7?(") 


Ai? (") 


/peak (erg/s cm^ "'^) 


F (erg/s cm^) 


CS/S-ratio 


Gunnz 


total shell 
star 






<3xl0~''' 


<2xl0"'^ 
2.1x10"** 


<lxlO~* 


KI 


shelll 


25: 


3: 


4.3x10"": 


5.3x10"'*: 






shell2 


36.5 


6: 


9.4x10"^**: 


3.6x10"'*: 






shell3 


43.2 


3.2 


2.9x10"^'' 


1.0x10"'^ 


2.4x10"* 




shell4 


46.1 


10.0 


3.8xl0"^'' 


1.2x10"^^ 


2.9x10"=^ 




total shell 








1.2x10"'^ 


2.9x10"* 




star 








4.2x10"^ 




NaD 


shelll 


25: 


8: 


7.6x10"'': 


1.9x10"''': 






shell2 


36.5 


6.4 


9.0x10"" 


3.5x10"'^ 






shell3 


43.4 


3.5 


1.2x10"^'' 


4.1x10"'^ 


2.3x10"* 




shell4 


46.1 


10.0 


2.0x10"" 


6.2x10"^* 


3.5x10"=^ 




total shell 








1.0x10"^^ 


5.6x10"* 




star 








1.8x10"^ 




b 


shell4 
star 


46.1 


10.0 


5.5x10"'** 


1.8x10"'^ 
4.8x10"'° 


3.7x10"* 



data of detached shells (Olofsson et al. 200(]|) . This gives 
a maximum CS/S-ratio of about 2x10""* and 4x10"^ for 
line scattering in the narrow KI ( only one of the doublet 
lines) and NaD (both doublet lines) filters, respectively 
(aeff~0.5 for the strong resonance lines, AAf=5nm). A 
comparison with the observed values. Tables || and |[ in- 
dicates that we are close to these ratios in both the KI 
and Na D filters for both sources (exceed them by at most 
a factor of four). For dust scattering the maximum CS/S- 
ratio is 1, since AAsc=AAf and aoff=l- 

A more strict limit is provided by the fact that the 
observed surface brightness in a filter cannot exceed the 
source function for scattering at the position of the shell, 
i.e., the mean intensity of the stellar radiation and the 
scattered radiation in the optically thick shell. This results 
in 

acff-F*,AAAsc 

where is the stellar flux density at the Earth averaged 
over the filter passband AAf, and 0sh is the angular size 
of the shell. 

Using the observed shell sizes the maximum surface 
brightnesses for line scattering are 5x10"^^ and 4x10"^^ 
erg/s cm^"^ for R Scl and 1x10"^^ and 1x10"^^ erg/s 
cm^ for U Ant in the KI and NaD filters, respectively 
[from now on, in the case of U Ant, we concentrate on the 



43" shell (shells) for which we also have information from 
the CO radio line observations]. The observed peak surface 
brightnesses are at most a factor of four above these val- 
ues. The maximum surface brightnesses for dust scatter- 
ing are much higher, 2x10"^^ and 1x10"^^ erg/s cm^ "^ 
for R Scl in the KI and NaD filters, respectively, and 
7xl0"i3, 3xl0"i3 and 7x10"" erg/s cm^"^ for U Ant 
in the KI, NaD, and b filters, respectively. 

We will proceed with a simple analysis of the observed 
scattered emission based on the assumption of optically 
thin scattering in a shell with no velocity gradient. The 
circumstellar scattered flux at the Earth, integrated over 
the total scattering wavelength range, in the optically thin 
regime is given by 

so that the expected CS/S-ratio is given by, 

Fac -Fsc,aAAsc AAsc 1 

-FT = -F: TTV = -XT' A U2 ^^scgsc«cff, (4) 

where Ngc is the number of scatterers, CTsc is the scattering 
cross section, and R^h is the shell radius [(^) is the opti- 
cally thin limit of Eq.(|l|)]. The assumption of optically 
thin scattering requires that the total scattering surface is 
smaller than the geometrical surface, i.e., 

NscCTsc ^ , 
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4-.1. Line scattering 

In the case of line scattering and assuming a Boltzmann 
distribution of the level populations, the number of scat- 
terers and the cross section are given by, 

El 

TlnfxMsh gie kT 



1 1 q'X'f 



(6) 
(7) 



where Mgh is the shell mass, is the mean molecular 
weight, Z is the partition function, gi is the statistical 
weight of the lower level, Ei is the excitation energy of 
the lower level, fx is the number abundance of species X 
with respect to H, rj^ is the fraction of X that remains in 
neutral form, and / is the oscillator strength of the line 
(for the other parameters we use common notations). 

It is now possible to use Eqs. (||), (|6|), and (^ to esti- 
mate whether the observed emission in the KI and NaD 
filters can be attributed to line scattering. Using the shell 
radii obtained from the observations, assuming a shell 
mass of 0.01 M0, which is a typical value derived from 
CO observations (Schoier & Olofsson 2001 ), and a solar 
abundance of potassium, and assuming a local velocity 
width of lkms~^, we find that in order to reproduce the 
observed CS/S-ratios in the KI filter only «6% of the K 
atoms in R Scl and ^4% in U Ant must remain in neutral 
form. The corresponding values assuming a solar abun- 
dance of sodium is «0.3% for R Scl and ^0.2% for U Ant. 

The transition between optically thin and thick line 
scattering can be obtained from Eqs. (||), (^, and (^. 
Saturation in the KI line is reached with 77n=0.02 for R Scl 
and ryn=0.04 for U Ant. For the NaD lines, saturation is 
reached already with 77n=6xlO^** for R Scl and 77n=10^'^ 
for U Ant. The values of rj^ required to obtain the observed 
CS/S-ratios in R Scl are higher than this by a factor of a 
few, and hence any KI and Na D line scattering is likely to 
be optically thick. In the case of U Ant, at least partially, 
optically thin KI and NaD line scattering is possible. We 
note here that these estimates are highly uncertain, and 
shall only be taken as a guide for the final interpretation. 

The ionization degree of the atoms in the shells is de- 
termined by stellar, as well as interstellar, UV photons. If 
we consider the stars to be blackbodies with a temperature 
of 2500 K and a luminosity of 5000 L©, and examine, as a 
first approximation, photons of 6 eV (this is in the region 
where the KI photoionization cross section has its maxi- 
mum, and also close to the ionization energy of Na), we 
find stellar photon flux densities at the shell radii of about 
2xl0~^ cm~^ s~^ Hz"^ in the case of R Scl and 10~^ cm~^ 
s""'^ Hz~^ for U Ant (note that this is very dependent on 
the assumed photon energy, and an accurate calculation 
requires an integration over the ionization cross section 
and the true stellar radiation field, which declines rapidly 



tons clearly dominates at the present time. The atoms, 
however, have been exposed to a much stronger stellar 
radiation field in the past when they were closer to the 
star. The stellar and interstellar photon fluxes were com- 
parable at a radius of RilO^^ cm. In other words, in R Scl 
and U Ant the stellar radiation dominated for about 10% 
and 5% of the expansion time, respectively. The ioniza- 
tion is counteracted by recombination, but it would lead 
too far in this connec tion t o analyse this in det ail (se e e.g. 
Glassgold & IIuggins [l98^ ). Mauron & Caux ( [l992[ ) have 
estimated that the ionization rate of KI due to the in- 
terstellar radiation field is about 5x10^^^ s~^. This gives 
an expected lifetime for the neutral species of Ri650yr if 
the interstellar radiation field dominates and there is no 
shielding. Taken into account that both the R Scl and 
U Ant shells are around (2-3)xl0'^yr old, the amount of 
K that remains in neutral form in these shells is «0.02. 
This number is in reasonable agreement with the above 
required values of rjn- We expect the ionization degree of 
Na to be comparable to that of K, and hence that the 
Na scattering is considerably optically thicker than the K 
scattering. 

Thus, we find that the emission in the KI and NaD 
filters can, in principle, be explained by scattering in the 
resonance lines of those potassium and sodium atoms that 
remain in the detached shells. Our simple calculations also 
suggest that the NaD lines can produce optically thick 
scattering (as seems to be the case for both R Scl and 
U Ant), while KI line scattering may be at least partly 
optically thin towards U Ant, and this would explain the 
different appearances in the U Ant KI and NaD filter 
images. 



4-2. Dust scattering 

In the case of dust scattering, and assuming single-sized, 
spherical grains, we have the following relations for the 
number of scatterers and the scattering cross section, 

A^ec = , (8) 

CTsc = QscTra^ , (9) 
where /dg is the dust-to-gas mass ratio, pg is the density 
of a dust grain, a is the grain radius, and Qsc is the grain 
scattering efficiency. If the condition for Rayleigh scatter- 
ing is fulfilled (i.e., a/A^0.16) we have 



12871^ 
3A4 



Re 



1 



(10) 



at these short wavelengths). Draine (197S) gives an inter- 
stellar photon flux density at 6 eV of about 2x10"'^ cm~^ 
s""'^ Hz"-'^. Therefore, ionization by interstellar UV pho- 



where n is the complex index of refraction of the dust 
particle relative to the surrounding medium. 

Using Eqs. (||), (||), and (||), typical values for the 
different parameters, (/dg=0.005, pg=2gcm~^, a=0.1/xm, 
Qsc<l), and A/sh=0.01MQ, we find for the shells of R Scl 
and U Ant that the area-ratios are <0.02 and <0.01, 
respectively, i.e., any dust scattering must be optically 
thin. Alternatively, the circumstellar medium is highly 
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clumped, which may indeed be the case (Olofsson et al. 
2000| ). 

We can now estimate whether it is possible to fit 
the CS/S-ratios in the KI, NaD, and b (for U Ant) fil- 
ters assuming only dust scattering. Using Eqs. (Q), 
and (^), and the fact that AAsc=AAf, we find that this 
is, in principle, possible provided that Qsc/a is about 
4xl0'^cm~"'^ and only weakly dependent on the wave- 
length. The latter requires a>0.1^m, and for such qrains 
(Qsc/a)S>4x 10^ cm~^. In addition, however, we have the 
requirement of optically thick emission for R Scl in the 
KI and NaD filters and U Ant in the NaD filter, which 
is obviously not consistent with the low CS/S-ratios and 
dust scattering, unless the circumstcllar medium is highly 
clumped. 

We also estimate whether it is possible to explain the 
CS/S-ratios for the component shell4, which has a very dif- 
ferent radial distribution than the bulk emission in the KI 
and NaD filters, by dust scattering. This is possible if we 
assume a shell mass of 0.01 M©, /dg=0.005, pg=2gcm~^, 
and Rayleigh scattering, i.e., we substitute (|l^) in @, 
with an average grain size of 0.03/xm and h=1.6. This 
grain size is consistent with Rayleigh scattering also in 
the b filter since a«0.05A. 

5. Discussion 

5.1. The nature of the scattered light 

No doubt, the circumstellar emission that we detect is due 
to scattering in atomic/molecular lines and/or by dust 
particles. We proceed now with a number of argument for 
and against different scattering agents. 

First, we note that the observed CS/S-ratios and sur- 
face brightnesses are remarkably close to the maximum al- 
lowed values for line scattering in the KI and NaD filters, 
and well below those allowed by dust scattering. We be- 
lieve that the uncertainties in the calibration of our data, 
and in some of the assumptions made (e.g., no velocity 
gradient, and the adopted local velocity width, shell mass, 
etc.), are such that a factor of a few can be accomodated 
with margin. 

Second, Gustafsson et al. (1997) have shown, using 



spectroscopic observations, that in the case of R Scl there 
exist KI resonance line scattered emission out to at least 
12" from the star, and that its brightness decreases as the 
third power of the angular distance from the star, i.e., they 
see no hint of a detached shell. However, beyond about 
12" the circumstellar emission is lost in the noise; their 
CAT/CES observations are clearly less sensitive than our 
3.6m observations due to the smaller telescope and the 
optical losses in the spectrometer. We find a brightness of 
about 2x10""'^^ erg s~^ cm~^ arcsec"^ for the circumstel- 
lar emission in the KI filter beyond wl2". This is a factor 
of five lower than the brightness obtained by Gustafsson et 
al. at an offset of wlO", and hence the data are not incon- 
sistent, although there is considerable uncertainty in our 



measured fluxes and to a lesser extent also in the fluxes 
reported by Gustafsson et al.. 

Finally, the U Ant images appear very different in the 
three filters. In particular, the difference between the KI 
and Na D filter images could be naturally explained as an 
effect of, at least partially, optically thin line scattering in 
the former, and essentially optically thick line scattering 
in the latter. Also the appearance of the R Scl KI and 
NaD filter images, i.e., the uniform surface brightnesses, 
can only be explained by optically thick scattering, and 
since the observed CS/S-ratios are much less than one, 
this requires a highly clumped medium in the case of dust 
scattering. 

On the other hand, the observed emission in the b filter 
has a very different radial distribution than those in the 
KI and NaD filter images, and an interpretation in terms 
of dust scattering appears more reasonable here. There are 
reasons for obtaining a positive result only in the b filter. 
For small particles the expected wavelength behaviour for 
the scattering process is such that we expect a higher con- 
trast between the scattered and direct stellar components 
in a blue filter than in a red one. For observations in a 
filter with a much redder central wavelength than that 
of the b filter (e.g., Gunnz) the reduction becomes more 
problematic, since the star is so bright that the circumstel- 
lar emission is swamped by the photospheric light even in 
the presence of a coronographic mask, and the scattering 
efficiency is also lower. On the other hand, in a filter cen- 
tred at even shorter wavelengths (e.g., v) the cool carbon 
stars have such a dramatic drop in their brightness (typ- 
ically, B-VRi3) that the number of stellar photons which 
can be scattered decreases drastically, pushing the circum- 
stellar emission down below the detection limit, even if 
the dust scattering efficiency increases further at shorter 
wavelengths. 

Although we lack a definite proof (e.g., no spectro- 
scopic information), we favour an interpretation where the 
main emission in the KI and NaD filters are due to res- 
onance line scattering. We further argue that weak fea- 
tures in the KI, NaD, and in particular b, filters are due 
to Rayleigh scattering by small grains. Similar observa- 
tions done in a filter containing no resonance lines and 
centred on the wavelength range between the NaD and 
KI lines (e.g.. Ha) will ultimately help us discern whether 
the emission is due to resonance line or dust scattering. 
Unlike the b and Gunnz filters, the stellar fiux and the 
dust scattering efficiency in such a filter lie in between 
their values at the NaD and KI wavelengths in which we 
have clearly detected circumstellar scattered light. Imag- 
ing through a very narrow filter which is tunable over the 
KI and/or NaD lines would also be an attractive observa- 
tional mode. Awaiting such data, our interpretation must 
be regarded as only tentative. 

This interpretation leads to different locations of the 
gas and the dust, in terms of the outer edges of shellS 
and shelly the difference is about 7" (i.e., -Rdusf-Rgas ~ 
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3x10"'^^ cm), which may at first seem surprising. However, 
this can be interpreted in terms of a drift of the dust par- 
ticles with respect to the gas, produced by the radiation 
pressure on the dust grains. Using a gas shell expansion 
velocity of 18.1 km s~^ as derived from CO observations 



(Olofsson et al. 1996), this drift velocity is estimated to 
be about (i?dust — Rgas)ve/Rga.s ~ 3kms^^. This is a per- 
fectly reasonable value for a mass loss rate approaching 
10~^ Mq yr"-'^ (Kwok 1975), which is in the range of mass 



loss rate estim ates for the detached CO shells (Schoier & 
Olofsson |200l| ). 



5. 2. The structure of the shells 

The images of R Scl and U Ant obtained in circumstellar 
scattered light confirm the results on detached circumstel- 
lar shells obtained previously from CO radio line data, but 
also reveal new interesting features. 

The CO data on R Scl show a marginally resolved de- 
tached shell. Modelling of the CO data, as well as HCN 
and CN radio line data, suggests a relatively thick shell, 
«12", with a mean radius of Ril2", i.e., an outer radius 
of ssl8" (Olofsson et al. 1996). This outer radius agrees 



reasonably with the relatively sharp outer 'disk' edges at 
R:!21" in our optical images, so that the radio line emis- 
sion and the scattered emission may be consistent with 
each other, but higher spatial resolution CO observations 
and more detailed modelling of the CO data are required 
to resolve this question. In particular, the inability to sep- 
arate the emission from the present mass loss envelope 
and the detached shell make the modelling of the radio 
lines very uncertain. The high optical depths of the cir- 
cumstellar scattered emission will make it difficult to ob- 
tain more information on the shell's structure from such 
data, although polarization observations may help in this 
respect. 

Additional information on the R Scl shell is obtained 
from the spectroscopic KI observations of Gustafsson et 



al. ( |1997| ). They found that their data, within 10" of the 
star, are consistent with a constant mass loss rate of about 
5x10~^Mq yr~^ if 100% of the potassium atoms remain 
neutral. This amounts to a mass of about SxlO^^'M© in- 
side 10", i.e., considerably lower than the estimated CO 
shell mass. This suggests that the shell has detached and 
that there is much less mass inside about 10", but the 
mass estimate of Gustafsson et al. can be higher if the 
amount of neutral potassium, for which there is presently 
no constraint, is lower. 

Using the expansion velocity of the detached CO shell, 
15.9kms^^ (Olofsson et al. 1996), and a distance of360pc 



derived from the period-luminosity relation for carbon 



stars (Groenewegen & Whitelock 1996), we estimate that 
the outer radius, 21", corresponds to an age of «2.3xl0'^ 
yr. At this time a drastic change in the mass loss prop- 
erties of R Scl must have occurred over a time scale of a 
few hundred years, as evidenced by the sharp cut-offs in 



the radial profiles of the scattered light. This conclusion 
is also valid in an interacting wind model, where the shell 
represents the shock zone between a fast and a slow stellar 
wind, but the time estimates may have to be redone. 

In the case of U Ant the CO data clearly reveal a 
detached shell whose radius of ks41" is relatively well de- 
termined, but where the width is unresolved at ssl5" res- 



olution (Olofsson et al. 1996). To our surprise, the optical 
images, in particular the one obtained in the NaD filter, 
suggest a multiple-shell structure, where the shell at «43" 
dominates in mass, and where there is at least one in- 
ner shell, at «37". The dominant CO emission apparently 
comes from the 43" shell. The fit of a shell brightness dis- 
tribution to this component in the KI filter radial bright- 
ness profile is not perfect, but nevertheless suggests that 
the shell is geometrically very thin, w3", which should be 
regarded as an upper limit since undulations in the shell 
structure, optical depth effects, and seeing will broaden 
the azimuthally averaged radial brightness profile. This 
width is comparable to those of the thin CO shells de- 
tected towards the carbon stars U Cam and TT Cyg 
(Lindqvist et al. 199£; Olofsson et al. 2000). We can also 



use the fit to the KI radial brightness distribution to esti- 
mate the masses of the inner shells relative to that of the 
outer shell. Taking into account that the number of pho- 
tons per unit area to be scattered decreases as the radius 
squared, we find that the masses of the two inner shells 
are only a few percent of the mass of the outer shell. The 
dynamic range of the CO maps is limited, and also con- 
sidering the poor spatial resolution, we can only conclude 
that any CO emission from the two inner shells must be 
at least a factor of five weaker than that of the 43" shell. 
Thus, the CO data do not exclude the existence of inner 
shells with the properties found here. The difference be- 
tween the morphology of the KI and Na D filter images can 
be explained by resonance line scattering and higher op- 
tical depths in the NaD lines. This could also explain the 
poor fit of geometrically thin shell brightness distributions 
to the NaD data. Optical depth effects tend to broaden 
the profile towards the star, while leaving the outer edge 
essentially unchanged. It is not obvious why we see shells 
inside the 43" shell if the Na scattering is optically thick. 
Two possible explanations are a clumpy medium or differ- 
ences in kinematics of the shells which will make the outer 
shell transparent to emission from the inner shells. 

Using the CO shell expansion velocity of 18.1 km s~^ 
(Olofsson et al. 1996| ), and a distance of 260 pc to U Ant 
derived from the Hipparcos parallax (Alksnis et al. 1998), 
the ages of the shelll, shell2, and shell3 components are 
estimated to be 1.7, 2.6, and 3.0x10^ yr. The time in- 
terval between the shells is consequently 400-900 yr. The 
morphology resembles to a large extent the multiple-shell 
structures (in the form of arcs, i.e., incomplete shells cen- 
tered on the star) seen towards the extreme carbon star 
IRC-t-10216 (time interval 200-800 yr; Mauron & Huggins 



1999i POOOD , and the post-AGB objects IRAS 17150-3224 
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(200-300 yr; K wok e t al. |1998| ), and CRL2688 (150-450 
yr; Sahai et al. 1998), and similar structures towards the 
M supergiant fi Cep where shells are seen in KI emis- 
sion (Mauron 1997). In these cases the shells represent 
density enhancements on a general density law. Sug- 
gestions for mechanisms for the formation of these shell 
structures range from modulated mass loss due to a com- 
panion (Harpaz et al. 1997[ ), shock waves in the CSE due 
to a close binary (Mastrodemos & Morris 1999), mass loss 



variations due to surface effects related to a magnetic ac- 
tivity cycle strengthened by a binary companion (Soker 
'Odj ), and pulsation induced effects internal to the star 
(Icke et al. 1992). It is not clear whether the putative 



multiple-shell structure seen towards U Ant is of the same 
type. An attractive alternative explanation is that there 
exist only two shells {shell2 and shellS), of which the outer 
one is the shock formed by a fast, short-term wind run- 
ning into a slower, long-term wind, and the inner shell is 
the reverse shock set up by this interaction. The shells are 
separated by 6" and over the lifetime of the outer shell 
this corresponds to a reverse velocity of about 3kms^^, 
which is not unreasonable. We will await more data (e.g., 
obtained in polarization mode) before further speculating 
on this tentative multiple-shell structure and its relevance 
for mass loss on the AGB. 
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